Geometric phase of an atom inside an adiabatic radio frequency potential 
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We investigate the geometric phase of an atom inside an adiabatic radio frequency (rf) potential 
created from a static magnetic field (B-field) and a time dependent rf field. The spatial motion 
of the atomic center of mass is shown to give rise to a geometric phase, or Berry's phase, to the 
adiabatically evolving atomic hyperfine spin along the local B-field. This phase is found to depend 
on both the static B-field along the semi-classical trajectory of the atomic center of mass and an 
"effective magnetic field" of the total B-field, including the oscillating rf field. Specific calculations 
are provided for several recent atom interferometry experiments and proposals utilizing adiabatic rf 
potentials. 
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I. INTRODUCTION 

Magnetic trapping is an important enabling technol- 
ogy for the active research field of neutral atomic quan- 
tum gases. A variety of trap potentials can be developed 
using magnetic (B-) fields with different spatial distri- 
butions and time variations. For instance, the widely 
used quadrupole trap and the loffe-Pritchard trap [l| are 
usually created with static B-fields, while the time aver- 
aged orbiting potential (TOP) 0] and time orbiting ring 
trap (TORT) ^3|, i^] are created using oscillating B-fields 
with frequencies larger than the effective trap frequen- 
cies. Atom chips have brought further developments 
to magnetic trap technology, as they can provide larger 
B-fields and gradients at reduced power-consumptions 
or electric currents using micro-fabricated coils. Today, 
magnetic trapping is a versatile tool used in many labo- 
ratories around the world for controlling atomic spatial 
motion in regions of different scales and geometric shapes, 
e.g., 3D or 2D traps, double well traps, and storage ring 
traps P,0,i,lllEi,l3- 

Recently, magnetic traps based on adiabatic microwave 
d, [1] and adiabatic radio frequency (rf) potentials 

]ARFP)Jio|JnL[il, [H, M M M mMMM HH, 

m, m, l24l l25l l26j have attracted considerable attention. 
An ARFP is typically created with the combination of 
a static B-field and an rf field. The idea for an ARFP 
has been around for some time [l^, [ll|, [l^l , and experi- 
mental demonstrations recently have been carried out for 
confiniiig both thermal and Bose condensed atoms 
[mil, lli, [13 • Further development of improved ARFP 
with atom chip technology likely will assist in practi- 
cal applications of atom interferometry. For instance, 
a double well potential was constructed recently using 
low order multi-poles capable of atomic beam splitting 
while maintaining tight spatial confinement |21]. Sev- 
eral interesting recent proposals outline the construc- 
tion of small storage rings with radii of the order 1/im 
[l8, 19, 20, 21, 22], which could become useful if imple- 
mented for atom Sagnac interferometry [27| setups on 
atom chips. 



When a neutral atom is confined in a magnetic poten- 
tial, its hyperfine spin is assumed to follow adiabatically 
the spatial variation of the B-field direction during its 
spatial translational motion. As a result of this adia- 
batic approximation, the center of mass motion for the 
atom experiences an induced gauge field [11], giving rise 
to a geometric phase (or Berry's phase) to the atomic 
internal spin state [2^. [3Qi]. The effect of this geometric 
phase is widely known, and is first addressed carefully 
in a meaningful way for atomic quantum gases by an 
explicit calculation of the resulting geometric phase in a 
static or a time averaged magnetic trap in Ref. [3l| . Sev- 
eral important consequences are predicted to occur for a 
magnetically trapped atomic condensate in a quadrupole 
trap, a loffe-Pritchard trap [IH [11], or a TORT based 
storage ring |33i]. To our knowledge, this geometric phase 
effect has not been investigated in any detail for an atom 
inside an ARFP. 

In a recent paper, we show that this geometric phase 
causes an effective Aharonov-Bohm-type [s^] phase shift 
in a magnetic storage ring based atom interferometer [33| . 
In addition, our studies imply that the spatial fiuctuation 
of the geometric phase can lead to a reduction of the 
visibility of the interference contrast. In view of this, we 
decided to carry out this study as reported here for the 
atomic geometric phase in an ARFP in order to shed light 
on the proposed high precision atom Sagnac interference 
experiment [27]. 

Analytical derivations for this study at some places 
become rather tedious and complicated. We therefore 
first will summarize our major results here for readers 
who may not be interested in the intricate details. We 
find that the geometric phase in an ARFP generally takes 
a more complicated form in comparison to the case of a 
static trap or a time averaged trap. In an ARFP, this 
phase factor is found to be determined by the trajectory 
of the time independent component of the trap field as 
well as an "effective B-field" that depends on the total 
B-field. In contrast to the earlier result found for a static 
trap or a time averaged trap [3l|,[3l], the final result turns 
out to be not expressible as a functional of the trajectory 
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for the direction of the total B-field in the parameter 
space. 

This paper is organized as follows. In sec. II, we gener- 
alize the semi-classical approach as outlined in Ref. [2l| 
for the operating principle of an ARFP to a form more 
convenient for discussing the geometric phase. Section 
III parallels that of sec. II by reformulating a full quan- 
tum theory for discussing the geometric phase inside an 
ARFP [24I. The explicit expression for the geometric 
phase inside an ARFP is given in a readily adaptable 
form for specifical calculations. In sec. IV, we discuss the 
effect of the geometric phase in several types of ARFP 
recently proposed for atomic splitters and storage rings 
[il [ia S, [Si [11 and beam splitters [H [2l|. Finally, 
concluding remarks are given in sec. V. 

II. A SEMI-CLASSICAL APPROACH 

In this section, we provide a semi-classical formula- 
tion for calculating the atomic geometric phase inside an 
ARFP. The semi-classical working principle for an ARFP 
is described in Ref. [2l|, although only for the special 
case when the atomic center of mass is assumed at a fixed 
location. In order to calculate the geometric phase, our 
formulation allows for the explicit consideration of atomic 
center of mass motion classically. In our approach, the 
geometric phase is obtained naturally, and the validity 
conditions for both the adiabatic and the rotating wave 
approximations are clearly shown for an ARFP. 

Inside an ARFP [U, the total B-field B{r,t) is the 
sum of a static field component Bs (r ) and an oscillatory 
rf field Bo {r, t) , which conveniently is expressed as 

Bo{f, t) = cos(wt) + ^ (f, t) cos{ujt + 77). (1) 

where f is the spatial position vector of the atom, uj is 
frequency of the rf field, and 77 is a relative phase factor. 

In this section, we will assume that the atomic spatial 
motion is pre-determined, i.e., f{t) is given (as a slowly 
varying function of time t). For weak B-fields, the system 
Hamiltonian is simply the linear Zemman interaction 

H[t)^gF^^BF■B[r{t),t], (2) 

where is the corresponding Lande g-factor and /is 
denotes the Bohr magneton. ?i = 1 is assumed. 

For a static or a time averaged magnetic trap, the 
Hamiltonian ([2]) varies slowly over time scales of the Lar- 
mor precession of the atomic spin in the total B-field. 



During the eifectively slow trapped motion, the atomic 
hyperfine spin is assumed to be fixed at the instanta- 
neous eigenstate of the Hamiltonian ([2]). The geometric 
phase then can be calculated straightforwardly from the 
variation of the B-field direction in the parameter space 

In an ARFP, the situation is more complicated. Al- 
though the variation of Bs[f(t)] remains much slower 
than the Larmor precession, the rf frequency oj usually 
is assumed to be nearly resonant with the precession fre- 
quency. Thus, the Hamiltonian ([2|) contains both fast 
and slow time varying components, making the direct 
calculation of the geometric phase a more involved task. 
In the following, we will proceed step by step, clarifying 
the various approximations adopted along the way. 

To understand the working principle for an ARFP, we 
first decompose the Hamiltonian H{t) into the fol- 
lowing form 

H{t) - Hsim] + i/+[r(t)]e— * + HAr{t)V^\ (3) 

where and H± are all slow varying functions of time 
and are given by 

HMt)] = 9F^^BF -BsW)], 

H+W)] = \9Ft^BF ■ {B^^\r{t)]+e-''B'f^W)]) . 
H.[r{t)] = HlW)]- (4) 

Hs is diagonal in the spin angular momentum ba- 
sis defined along the local direction of the static B- 
field Bs[r{t)]. The eigenstate takes the familiar form 
\'mF[r{t)]) quantized along the direction of Bs[f[t)\, 
with the eigenvalue to^ |i3s[r(i)]| for Bs[r{t)] ■ F and 
rap £ [—F,F], in analogy with the usual case of the z- 
quantized representation result of Fz\mp)z — mp\mp) z- 

Next we introduce a unitary transformation 

F 

U{t)= \mp).s{mp[mW^'''''\ (5) 

mF— — F 

with K = sign(5i?) for the rotating wave approxima- 
tion. The quantum state in the interaction picture 
\^{t))i = U{t)\^{t)) defined by U{t) is governed by the 
Schroedinger equation idt\^{t))i = Hi{t)\^{t))i, with 
the Hamiltonian in the interaction picture given by 



Hi{t) = UHU^ +i{dtU)U'^ 

F F , 

= J2 ™«;A[r(t)]|m),,(m|-z ^ |m), ,(™[f(i)]|||n[f(t)]), ,(n|e^(™-")«-* 

m— — F m,n— — F 
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{hli\m]H. Am~l\ + /iL")[r(0]|m), ,(m - lle^'""* + h.c.) 



m=-F+l 



m z z{m e 



h.c) 
I 



(6) 



where the time dependent parameters are defined as 

A[r(i)] = ^XB\gFB[rm-^, 
h^^\r{t)] ^ ,(™[f(t)]|ff±[rW]|(m-l)[r(i)]>„ an^) 
h„r[r{t)] - s{m[r{mH±W)]Mrm)s- 

The above result is obtained easily if we note that 
the matrix element s{m[r{t)\H±{t)\m'[r{t)) s is non-zero 
only when m — m' = Q, ±1. So far, we have always as- 
sumed that |m[r(t)])s is a single valued function of the 
atomic position r. A careful examination shows that the 
eigenstate |m[r(t)])s cannot be determined uniquely be- 
cause of the presence of the U{1) gauge freedom for se- 
lecting a local phase factor exp{i0[r (<)]}, which conse- 
quently affects the resulting expressions for hm \t) and 
,{m[r{t)]\d/dt\m'[f{t)])s. 

The rotating wave approximation neglects of the 
oscillating terms proportional to e*""^* (m 7^ 0) in the 
Hamiltonian Hj ([6]). The error for this approximation 
is estimated easily from a time dependent perturbation 
calculation. The sufhcient condition for its validity 
requires that all factors such as dt'hm{t') exp[iKUJt'], 

J*dt'hi;;\t')^m,m-i{t')eMM'2^ + ^)t'], and 
Jl^dt'{m[r{t')]\d/dt'\n[r{t')])^,nn{t')ciip[i{rn - n)K{LO + 
A)t'] are negligible, where 



Cmn(i) = exp 



dt' ,{m[f{t')]\ — \m[rit')]) 



exp 



dt' s{n[f{t')]\-\n[rit')]) 



(8) 



Thus, the gauge independent factors /im^^m,m-i, 
{'ms\d/dt\ns)^mn , h„n and A should all vary slowly with 
time and with the modulus of their amplitudes much less 
than Lo. 

The effective Hamiltonian in the interaction picture 
under the rotating wave approximation then becomes 

F ^ 

~i E \^)^s{m[r{t)]\-\m[f{msz{m\,{9) 



m— — F 



dt' 



where the first term resembles a coupling between the 
atomic spin and an "effective B- field" B'^^{r), whose 
components in real space are given by 



Re 



2,(m(f)|H±(r)|(m-l)(r)), 
^iB9F^/{F + m){F -m-l) 



^cff 



B 



y (r) 



Tni 



2s{m{f)\H±{r)\{m-l)if))s 



fiB9F y/{F + m){F -m-1) 



Bs{^ 



and 



(10) 



^^B\gF\ 

Clearly, the x- and y-components of the effective field 
B°^{f) depend on the explicit form of the eigenstate 
|m(r))s. In fact, it easily can be seen that different 
choices of the local phase factor for the \m{r))s actually 
lead to different values of B^^{r) related to each other 
through r-dependent rotations in the x-y plane. 

In practice, the eigenstate \n{r))s and the effective field 
B'^^ can sometimes be constructed more simply, as in 
Ref. [2l[. For any spatial position r, we first choose a 
rotation R[rh{r) , xi'f^] along the axis m(r) with an angle 
x(r) that satisfies R[n{f),xi'f^]Bs{r) = \Bs{r)\ez. It is 
then easy to show that the eigenstate |n[r(i)]),; can be 
chosen as 



|n(f))s=exp iF ■ 7fi{f)x{r) \n)z. 



(11) 



Unfortunately, the choice for R is not unique in a given 
static field Bs{r), an analogous result to the U{1) gauge 
freedom for the the egienstate \n{r))s- Corresponding 
to the choice ([TT]) given above for |n(r))s, the unitary 
transformation U defined in ([5]) would become 



U{t) — exTp{—iFzLut) ■ exp —iF ■ m^rjxir) 



(12) 



and the transverse components of the "effective B-field" 
given by B^l{r) = B,^y{f)/2 |21j with 



Bir) 



R[m{r1,XmBlf{r^ 
+R[e,,~Kr,]R[mir),xir^]Bi'\r^. (13) 

In earlier discussions of an ARFP [U [S^l , the atomic 
internal state is assumed uniformly to remain adiabati- 
cally in a certain eigenstate of the first term of Hj^^ (t) . 
To fully appreciate this adiabatic approximation and to 
calculate the geometric phase, we expand |^'(t))/ into 
the instantaneous eigenstate basis |n[r(t)])oft quantized 
along the direction of the effective B-field J5°* accord- 
ing to |*(i)>/ = EnCn{t)\n[f{t)])cS. The first term of 
H^g{t) is simply the effective Zemman interaction be- 
tween the atomic hyperfine spin and the effective B-field. 
The corresponding Schroedinger equation for the Hamil- 
tonian H^g (t) of ([9]) then becomes 



iJ^Gn{t) = [er^(0 +^n„(t)]a.(i) + E ''nm{t)Cmit), 
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with 



I 



-*eff(p[r(i)]|^k[r(i)])eff. 



(14) 



Under the adiabatic approximation, the atomic inter- 
nal state remains in a given eigenstate |n[r(t)])eff with 
transitions to states |m[r(t)])eff (m ^ n) being neghgi- 
bly smaU. Thus, the transition probabihty, as estimated 
from the first order perturbation theory, 



should be much less than one. As before, we find the 
sufficient condition for the validity of the adiabatic ap- 
proximation is given by 



\Vynn{t)\ 



\er{t')-^^P\t')\ 



« 1, 



(15) 



provided that i^mn(i) exp[i [i^„„(t") - Vrnm{t")]dt'], 
which is independent of the local phase factor for 
|n[r(i)])eff and |n[r(i)])s, remains a slowly varying func- 
tion of time. 

A straight forward calculation from the effective 
Hamiltonian ([9]) then gives the general expression for the 
geometric phase in an ARFP 



ln{t) 



ynn[t')dt' 



Y.\Mr{t') 



W{f)]\-i^W{t')])sdt' 



(16) 



and -/i'\t) = -I dt' ^s{n[r{t')\d/dt'\n[r{t'))^g. During 
the adiabatic motion in a given internal state, the time 
evolution of the coefficient C„ (t) takes the form 

C„(i) = C„(0)e-^^'^^"'(*')'^*'e-^T"(*). (17) 

Equation (flB)) is the central result of this work. The 
geometric phase of an atom inside an ARFP is shown to 
contain two parts. The second part, ^n\t), is clearly 

due to the interaction term fisgpF • in H^^^ 
with its value determined by the trajectory of the direc- 
tion for the "effective B-field" B'^^ . The first part arises 
from the second term of hI,^ ([9]). It is determined by 
the trajectories of both the static field Bg and the ef- 
fective B-field B'^^ . The expression for 7„ in an ARFP 



is complicated because the internal quantum state in an 
ARFP is assumed to be adiabatically kept in an eigen- 
state of iJugpF ■ B"^^ , rather than an eigenstate of the 
total interaction Hamiltonian -ff^ff (^) ■ 

In section IV, we will perform explicit calculations for 
several examples of ARFP proposed for various applica- 
tions: e.g., as atomic storage rings or atomic beam split- 
ters. Most often we find that only the first part of Eq. 
(I16p contributes a non-zero value to the geometric phase. 

Before proceeding to the next section for a quantal 
treatment of the geometric phase, we find the time evo- 
lution of the atomic spin state in the Schroedinger picture 



\^it)) = 5]CK0).(m|/[r(t)])eff 



ml 



e "Jo 



«/o f-\[t')dt' -i-ii(t) -imut 



m[r{t)])s, (18) 



obtained directly from \^{t)) = UHt)\^{t))i after the 
applications of the rotating wave and adiabatic approxi- 
mations. When the atom is prepared initially in a specific 
adiabatic state |n[r(i)])eff of the interaction picture, we 
arrive at the simple case of C;(0) = 5;„. 



III. A QUANTUM MECHANICAL 
TREATMENT 



In the previous section, we provided the result for the 
geometric phase 7n(i) in an ARFP based on a semi- 
classical approach, where the atomic center of mass mo- 
tion is described classically. A clear physical picture ex- 
ists in this case for the appearance of the geometric phase 
in a certain parameter space. The validity conditions for 
the rotating wave and the adiabatic approximations as 
obtained above are all formulated in terms of gauge inde- 
pendent forms. However, if the influence of the geomet- 
ric phase on the atomic spatial motion is to be included, 
e.g., as in the Aharonov-Bohm-typc, phase shift, inter- 
ference arrangement in an atomic Sagnac interferometer 
discussed earlier [s^, we would need an improved de- 
scription where both the atomic spin and its center of 
mass motion are treated quantum mechanically. 

In a full quantum treatment of the atomic motion, the 
quantum state of an atom can be expressed as \'^(t)) = 
^ 0;(r, i)!/)^, where 4>i{r,t) is the atomic spatial wave 
function for the internal state \l)z of Fz. The state \^{t)) 
then satisfies the Schroedinger equation governed by the 
Hamiltonian 



n^ — +gF^lBF■B{r,t), 



(19) 



with P being the kinetic momentum and M the atomic 
mass. 

The rotating wave and adiabatic approximations can 
be introduced now by defining the interaction picture 
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with the unitary transformation 

F 



U{t) 



F 



.n=-F 



(20) 



The state in the interaction picture = 
now is governed by the Schroedinger equation with the 
Hamiltonian Ticg — lAT-iW . Under the rotating wave 
and adiabatic approximations, we neglect transitions be- 
tween states \m)z and \n)z (m 7^ n) as well as the rapidly 
oscillating terms. We then obtain 

Hoff w \n):, z{n\'Hcs\n)z z{n\ 

n 

« E^a^l">-H, (21) 



(n) 

where the adiabatic Hamiltonian H)^^ for the n-th adia- 
batic branch is defined as 



ad 2M 



+ 4")(r) 



(22) 



with the effective gauge potential 

Mr) = -^J2\Mrmfsmmir))s 

I 

-ieff(n(r)|V|n(f))eff. (23) 

In this form, it is well known that the geometric phase 
7„ can be expressed as the integral of the gauge poten- 
tial An along the spatial trajectory for the atomic center 
of mass in an ARFP, i.e., one would expect generally 
that 7,1 = / An ■ dr. Similar to the result of the semi- 
classical approach, the gauge potential A„(r) can be ex- 
pressed as the sum of two parts. The first part in Eq. 
(|23| is the weighted sum of the atomic gauge potential 
— i s(Z(r)|V|Z(r))s from the static field Bs, while the sec- 
ond term is the atomic gauge potential from the "effective 
B-field" B""^. 

A full quantum treatment for atomic motion in an 
ARFP has been attempted earlier In fact, many of 
our formulations are identical to the results of Ref. [22] . 
For instance, it is easy to show that the unitary trans- 
formations Us, Ur, and Up in [12] are related directly 
to ours as UpUlfUl = U. The only difference concerns 
the gauge potential An that was neglected in Ref. [2^ . 
Thus, they did not give the expression for the gauge po- 
tential, and the result for the geometric phase was not 
obtained either . Our study shows that the neglect of 
the adiabatic gauge potential potentially can give rise to 
a final result, dependent on the choice of the local phase 
factors for the internal eigenstate. 



IV. GEOMETRIC PHASES IN ARFP BASED 
APPLICATIONS 

In the above two sections, we obtain the expression for 
the atomic geometric phase in an ARFP. This section is 
devoted to the calculations of the geometric phases for 
several proposed applications of ARFP, such as storage 
rings or beam splitters for neutral atoms [l^, [l^, HO) El, 

Before presenting our results for the more specific 
cases, we provide some general discussions of the geo- 
metric phases in several ARFP based storage rings. As 
was pointed out earlier, the geometric phase 7„ is given 
by the line integral of the gauge potential An along the 
trajectory for the atomic center of mass motion. For a 
closed path in the storage ring at a fixed p — pc and 
z = Zc, this can be further reduced to 



In = q 



(24) 



where the integer q is the winding number of the path 

and A^f ^ is the component of An along the azimuthal 
direction of the familiar cylindrical coordinate system 
(p, (j>, z). Without loss of generality, we take g = 1 in this 
paper. For the storage rings proposed in Refs. [H, [H, 

HH, 112] ' tl^s gauge potentials A^\p, 0, z) are actually 
independent of the angle (j). Therefore, the geometric 
phase is simply given by 



in 



27rp,4«(pe,Ze), 



(25) 



given out in explicit forms for different storage ring 
schemes [H, [H [l^, [H [H . 

In reality, because of thermal motion or when the 
atomic transverse motional state is considered, the center 
of mass for an atom can deviate from (pc, Zc) even for a 
closed trajectory. This uncertainty in the exact shape of 
the closed trajectory gives rise to a fluctuating geometric 
phase and is usually difficult to study. Assuming a simple 
closed path at fixed p and z, we have found previously 
that the subsequently fluctuations could decrease the vis- 
ibility of the interference pattern , 33;] . Quantum mechan- 
ically, such destructive interference can be explained as 
resulting from entanglement between the freedoms for cj) 
and {p, z) because of the dependence of the gauge po- 
tential Af^ on p and z. Therefore, it is important to 
investigate this dependence near the trap center. 

For simplicity, our discussions below will focus on the 
closed loops where p and z are ^-independent constants. 
In this case, the geometric phase can be expressed as 
7„(p, z) = 2TTpAif\p, z). We will show numerically the 
distributions for 7n(p, z) obtained this way near the cen- 
tral region of {pc,Zc)- If needed, a more rigorous ap- 
proach can be developed to investigate the fluctuations 
of the resulting geometric phase from the gauge potential 
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A. The storage ring proposals of Refs. [l8l . Il9l . |2C 



This subsection is devoted to a detailed calculation of 
the geometric phases for the ARFP storage ring propos- 
als of Refs. [H, [H, [13] ■ We will derive the analyti- 
cal expressions for the azimuthal component Aif' of the 
gauge potential that arises in both cases from cylindri- 
cally symmetric static B-field and rf fields. Because of 
the cylindrical symmetry, the angle (3s{p, z) between the 
local static B-field and the z-axis is required to be ana- 
lytical in the region near the storage ring. Therefore, the 
eigenstate \n{r))s of F ■ Bs can be chosen as 

\n{f))s = exp{-i[F ■ e^l3s{p., z) + (26) 

Consequently, B°^{f) is also cylindrically symmetric, 
which leads to the eigenstate |n(r))cff of F ■ _B°^ as 

\n{r))otf - cxp{-^[F • nf{r)f3,sip, z) + 7i(^]}|n)^, (27) 

with the unit vector n'^ (r) in the x-y plane orthogonal to 
B'^^{r) and /3off(p, z) denoting the angle between B°^{r) 
and the z-axis. We note that the unit vector field ii'^{r) 
also possesses cylindrical symmetry, i.e., remains invari- 
ant under rotation around the z-axis. The expressions of 
(PS)) and ([?f|) allow us to obtain the simple expression of 
the gauge potential 

A^^\p,z) = --cos/3cff(p,2:)cos/3s(p,z), (28) 
P 

after straightforward calculations. 

In the scheme of Ref. [l^ , the static B-field is a "ring- 
shaped quadrupole field" that vanishes along a circle of 
a radius po in the x-y plane. Near p — po, the B- field is 
given approximately by 



Bsif) ^ B'{p- pa)ep- B'ze, 



(29) 



like a quadrupole field, while the rf-field takes a compli- 
cated form 

Bo{r,t) = cos(cji) -I- cos(cjt + (^)^ Cp 

— — ^ sin(a;t) + — ^ s\n{Lot + if>)^ Bz, (30) 

with constants a and b independent of f. 

From the expression of (pS)) for the eigenstate \n{r))s, 
the "effective B-field" becomes 



i?^«(r) - SV(p-Po)2 + z2-ro]e. 

6 ,^ . a 
—r= cos[0 + ip) -\ — j= cos t/ I e„ 

-^sin(6i + (^) - -^sin6l ) 60, (31) 




FIG. 1: (Color online) A cross-sectional view for the storage 
ring of Ref. The static field is zero in the ring at the 

fixed radius po- The addition of rf-fields creates an ARFP 
centered at a ring through {pc,Zc). The distance from the 
trap center to the ring with radius po in the plane z — is ro- 



where tq and 9 are given by 



cos6{p, z) 
sin0(p, z) 



IfiBgpB'l 

P- Po 

ViP-Po)^ + z^ 
z 

ViP-Po^ + z^ 



(32) 



In an ARFP, as discussed here, the trap center 
at {pc Zc) is determined by minimizing both the z- 
component and the transverse component of i?''^. With- 
out loss of generality, we will assume a,b > 0. Then, 
(pc Zc) is found to satisfy 



V {Pc - Pof + zl = ro. 



(33) 



i.e., the trap center lies on the surface of the "resonance 
toroid" at p = Po with a radius tq as shown in Fig. [1] 
The relative angle of the trap center with respect to the 
center of the toroid cross-section is given by — <p/2. On 
this "resonance toroid," the rf-field is resonant with the 
static field, i.e., B°^ vanishes. As a result, the "effective 
B-ficld" lies again in the x-y plane on the "resonance 
toroid," which gives cos/3eff (Pc, ^c) = and leads to the 
result A^\pc^ Zc) = 7„ = as shown in the t rap center 
for the storage ring considered before in Ref. |18l |. 

From the expression ([29| of the static field and the def- 
inition of the angle 9(p, z), we find a simple relationship 
Ps{Pi z) = Tr/2 -\- 9{p, z), with which the gauge potential 

a'^*\p,4>) in (Ull) can be further simplified as 
A':^\p,z) = -cosPes{p,z)sm9{p,z) 
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« -cosPcs{p,z)sme{pc,Zc), (34) 
P 

near the trap center. Thus, the spatial fluctuation of 
the gauge potential A^t"^ (p, z) in the region around the 
trap center is closely related to the angle 0{pc, Zc) of the 
trap center, or the parameter ip of the oscillating field 
Bo- When = 0, the atom is trapped in the region 
with « or TT, where the fluctuation of A^n^ {p, z) is 
suppressed significantly due to the small value of sin 6'. 
On the other hand, if the angle ip is set to tt with the 
trap center located in the region with 9 « ±7r/2, the 
fiuctuation of the gauge potential becomes amplified. 

In Fig. [21 we illustrate numerical results for the dis- 
tribution of the geometric phase 71 (p, 2;) = 27rpA^(p, z) 
in the region near the trap center &t Lp = 0, 7r/2,7r. We 
see clearly decreased fiuctuations of 71 when the absolute 
value of sin0(/9c, Zc) — — sin((/?/2) is decreased. 

Next we turn to the storage ring of Ref. [lH con- 
structed from a quadrupole static B-field Bs{r) = 
B'{x,y,—2z) and an r-independent rf field Bq = 
Bii cos{iut)ez along the z direction. The resulting ARFP 
provides a 2D ring shaped trap in the x-y plane. In ad- 
dition, a ID optical potential along the z direction is em- 

ted to confine atoms in the transverse plane at 2; = 
The "effective B-field" takes the form 

(f) = B'ip - po)e, - ^B,tep, (35) 

in the plane at z = 0, with po — uj/\pBgFB'\. Because 
the strength of B'^^ is near minimum at the ring p = po, 
the trap center for this storage ring is located at pc = 
Pq and Zc = 0. At the trap center, the "effective B- 
field" is along the direction of Cp. Thus, according to Eq. 

(|28|1 . the geometric phase 'jn at the trap center again 
vanishes. 

In Fig. [3l we show the distribution of the geometric 
phase 71 in the region near the trap center for _Brf = 
0.05|B'|po and B,i = 0.15\B'\po. We see that the fluctu- 
ation is relatively small when the strength of the rf-field 
is large. This can be explained by Eq. ([28)) , which shows 
that A^f^ is proportional to cos/3eff and can be approx- 
imated as 2B^^ /B^i near the trap center. When Bj-f is 
large, the gauge potential becomes a relatively slow vary- 
ing function of p and z. In this case, the presence of a 
ID optical potential allows for the possibility of tuning 
the trap center position to a nonzero value of z, with the 
storage ring remaining in the x-y plane. Then cos f3s is 
assumed to a nonzero value, leading to increased fluctu- 
ations for the geometric phase. 

Finally, we discuss the geometric phase in the "time 
averaged" ARFP storage ring proposed in Ref. [2^ . Un- 
like previously considered ARFP based storage rings, the 
time dependence now exists in both the "static B-field" 
and the frequency of the rf field given by 

Bs{r, t) = B' pCp - 2B'zez + Bm sin(cjmi)ez. 



(a) 




FIG. 2: (Color online) The distribution of the geometric phase 
71 near the trap center (pc,2c) of the storage ring proposed 
in Ref. ilS] at (a) = Q, (b) ip — it/2, and (c) <p = tt, clearly 
displaying the sm{(p/2) dependence. 



Bo{t) = B^{Sm[uj{t)t]ez, 
uj{t) = ^^0^1 + [B^/B'poY sin^{Lo^t) . (36) 

The frequency Um is assumed to be much smaller than cjq 
but much larger than the trap frequency. The radius po 
is now defined as po = wo/|p_b(7_f-B'|, and the "effective 
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(a) 



0.2 T 




FIG. 3: (Color online) The geometric phase 71 for the storage 
ring of Ref. ^ with (a) Brf = 0.15B' po and (b) Brf = 
O.OSB'po- 

B-field" takes the form 

B'«{r, t) = A(f, t)e, - i^'^ B,,e^ . (37) 
|2S,(f,t)| 

The operating; principle for the time averaged storage 
ring of Ref. [20] is similar to the well-known TOP [4I 
and TORT traps [1, 01 • The effective trap potential ex- 
perienced by the atom is proportional to the time aver- 
aged value of the "effective B-field" J^'''"'"' \B''^{r,t)\dt. 
When Brf and Bm are much smaller than B' p^, the 
center of the storage ring is located approximately at 
Pc ~ Poi — Q. Using the earlier result [33|, we find that 
in the time averaged storage ring, the effective gauge po- 
tential a'^'' (p, z) is reduced simply to the time averaged 
instantaneous gauge potential 

A'^^^ (P, = ^ r'^^^ (p, z, t)dt, (38) 

with A\f^ [p, z, t) given in (|28p . The geometric phase then 

is given approximately by 7n(p, z) = 2'KA\f^ (p, z). In this 
case, we find that the geometric phase always vanishes at 



(a) 




FIG. 4: (Color online) The geometric phase 71 for the storage 
ring of Ref. 20] at (a) B^i = O.SB'po and (b) Brf = O.lB'po- 
B,„ = 0.05B'po. 

the trap center (pc, ^c)- Figure [4] illustrates the distribu- 
tion of the geometric phase in the region near the trap 
center for two different values of the rf-field amplitude 
B-cf. Similar to the storage ring of Ref. [l^, the fluctua- 
tion of the geometric phase is suppressed in this case for 
large i?rf- 

B. The storage ring proposals of Refs. |2li. |22|| 

Next we consider the ARFP based storage ring pro- 
posed in Refs. [2l|, [S^. In this case, the static B-field 
is that of a loffe-Pritchard trap on an atom chip. In the 
Cartesian coordinate {x^y, z), it takes the form 

Bs^B'xe,-B'yey + B'Le,, (39) 

where B' is the B-field gradient and the bias field along 
the z-direction is denoted as B'L. The amplitudes i?^^"* 
and B^i^iz) of the rf field are B^^^ = [B,i{z) / y/2]e^ and 
4^f^ = [B,f{z)/V2]ey with 

B,,{z) = B';f +B"z^. (40) 
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In the schemes of Ref. [2l|, |22| considered earher, the 
phase r] of the rf field is assumed to be ktt/2. The x- 
and y-components of the "effective B-field" B°^{r) then 
become 



B 



cff / 



according to Eq. ([T 



0, (41) 

Then the strength of the "effective 
B-field" has its minimum along a circle with a non- 
zero radius pc, provided a positive detuning A exists at 
the origin (0,0,0) '21^ 21]- The "effective B-field" 5°* is 
easily shown to lie in the x-z plane along the trap bottom 
mapped out by the atomic center of mass motion. This 
gives rise to a vanishing ■ With a proper choice for 

the local phase of \n{r))cs, the gauge potential Alf'^ takes 
the form 



■cos/3off(p,2;)(l-cos/3s(p, z)). (42) 



P 

Figure [5] displays the geometric phase along a closed 
path for a spin-1 atom as a function of pc for the ARFP 
storage ring proposed in Rcfs. 21, 22]. The parameter A 
is defined as 



A = \/2- 



A[f= 0] 



(43) 



To assure the validity of the rotating wave approx- 
imation, we find that the maximal values of A[f — 

0]/(|5_f|a*b) and B^f^/V2 must be restricted to the re- 
gion of A G [0,0.15]. 

As shown in Fig. [l] the geometric phase remains much 
smaller than 27r in this situation. This fact can be appre- 
ciated easily if we look at the distribution of the "effec- 
tive B-field" B'^^. According to Eq. (gT]), the component 
B^^ has a nonzero minimal value i3rf/(2-\/2), while \B^^\ 
can become arbitrarily small, although not necessarily 
zero in general. Therefore, at the trap center where ji?''*^] 
is a minimum, the value of cos/3cff — B^^ /\B''^\ can be- 
come very small, leading to small geometric phases. Yet, 
despite the relatively small geometric phase found here, 
our result remains important because it could represent 
a systematic error if not properly included in a Sagnac 
interference experiment. 

In Fig. [SI we show the spatial distribution of the ge- 
ometric phase 7i around the trap center with A = 1/3 
and A = 3. The fluctuation is found to be relatively 
small when A is small or when the rf-field amplitude Brf 
is large. 

Although not discussed in Refs. [HI, [12], a ring shaped 
trap also can be realized if we take rj = —kit/ 2. The 
"effective B-field" B^^ still lies in the x-y plane 

B^z) 



B 



2V2 

Brijz) 

2V^ 



cos{2(f)){l ~ cos/3s(/9, z)), 
sin(2(/>)(l - cos/3s(p, z)). 



(44) 



r 

-0.02 
-0.04 
1 -0.06 
-0.08 
-0.1 



?i=3 



c) 



?i=1/3 \ 



0.1 0.2 0.3 0.4 0.5 0.6 

P /L 



FIG. 5: (Color online) The geometric phase 71 is plotted 
against the radius pc for the ARFP storage ring of Ref. 
[23, Hi] with 7? = KTi/2 at A = 3, A = 1, and A = 1/3. 
To assure the validity of the rotating wave approximation, in 
the solid lines, the maximal value of A[r = Q\/\gF p-bB' po\ 
or po) are restricted to be smaller than 0.15. The 

extending dashed line is beyond the rotating wave approxi- 
mation for A = 1/3 and B^f/{^/2B'po) G [0.15, 0.3]. 



clearly giving rise to a non-zero solid angle with respect to 
a closed path along the storage ring. Therefore, the term 
7^^'' is non-zero in this case. We choose the eigenstates 

\n{r))s and ln(f))off as 

\n{r))s = ex];>[-iF ■ ns(r)(3s{p,z)]\n)^, 
\n{r)),s - exp[-iF-ni*f(f)/3eff(p,z)]ln),. (45) 

with the unit vector n^(r) in the x-y plane orthogonal 

to Bs{r). In this case, the gauge potential A^'n'^ can be 
expressed as 



A^^\p. 



■ cos/?eff(p, z)[l + cos/3s(p, z)]. (46) 



In Figure [71 we show the fluctuation of the geometric 
phase 71 for a closed path with a new parameter 



A' = 6%/2 



A[f = 0] 



(47) 



equal to 3 and 1/3. The fluctuation for 71 is found to 
be much larger than the case of 77 = K7r/2, which can be 
explained by the transverse components i? ° ^ of the "ef- 
fective B-field." Because cos Ps is always close to unity. 
In the case of 77 = — K7r/2, B^^y can take only small pos- 
itive values. Therefore, at the minimum of the ARFP 
trap p = po of \B°^\, both Bf and B^^y have to be 
close to zero. In this case the value for cos /3eff becomes 
a rapidly changing function of p in the region near pc- 

Our above calculations have obtained analytical ex- 
pressions of the geometric phases in an ARFP based stor- 
age ring for rj = ±K7r/2. We have further investigated the 
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FIG. 6: (Color online) The spatial distribution of the geo- 
metric phase 7i for the storage ring of Refs. [2ll . [2^ at (a) 
A = 1/3 and (b) A = 3. 77 = k-k 12. Bf^ = O.OSB'L and 
B" = 10""B'/L are assumed. 



fluctuations of the geometric phase for the two cases of 
T] — -iiK,'n/2. It seems one benefits from implementing 
a Sagnac interferometer in the discussed ARFP storage 
ring with r\ — nix jl and operating at a relatively large A. 

Before proceeding onto the concluding section, we will 
discuss the geometric phase in an ARFP based beam 
splitter created via a double potential [3, HH ■ In such an 
implementation, the static field Bs is created from a loffe- 
Pritchard trap, while the oscillating rf field components 
are = i?rf[-z]ex and Bf^ — 0. By spatially tuning 
the amplitude of i?if from zero to a significant value, 
in the x-y plane, an ARFP can be tuned from a single 
well centered near the origin to a double well with two 
minimal points at the point with nonzero radius po a-nd 
(/) = 0, TT. Therefore, a Y-shaped atom beam splitter 
can be accomplished when the i?rf [^] initially is increased 
along the z-axis to a large value, and then decreased to 
zero. In such an arrangement, the atom beam moving 
along the z direction can be separated into two beams 
that move along the z-axis at = 0,7r for a while, and 
then can be recombined again into a single beam. 



(a) 




FIG. 7: (Color online) The spatial distribution of the geo- 
metric phase 71 for the storage ring of Refs. [2ll . [2^ at (a) 
A' = 3 and (b)A' = 1/3. t] = -/t7r/2. Bf^' = O.OSB'L and 
B" = 10~^^B'/L are assumed. 



In the atom interferometer considered above, both the 
static field B^ and the "effective B-field" B are limited 
to the x-z plane. Therefore, for motion along the closed 
path of the trap bottom, the solid angle enclosed by the 
trajectory of B"^^ is zero. Thus, the geometric phase in 
p6)) can be expressed as 



7«(t) = -^ / V|cff(n(r)|0.|^(^(^1|V|/(r)),•^/d^'. (48) 

Jo ; 



We can show that the product 
cff(n(r)|^)zp s(Z(r)|V|^(r))s is a function of pc and 
is independent of z. Thus, the geometric phase can be 
expressed as an integral of this function with respect to 
Pc, from zero to a large value and then back to zero. 
Therefore, the value of the geometric phase would be 
zero in the end. 
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V. CONCLUSION 

In this study, we develop theoretical formalisms for 
the calculation of the atomic geometric phase inside an 
ARFP. We show that, due to the complexity of the 
ARFP, the geometric phase depends on the spatial vari- 
ation of both the static field and an "effective B-field" 
_B°*^. We provide general expressions for the geometric 
phase and the corresponding adiabatic gauge potential in 
Eqs. p6|) and ((23)) . respectively. 

To shed light on actual applications of the atomic ge- 
ometric phase, we investigate the distribution of atomic 
geometric phases for several proposed or ongoing exper- 
iments with ARFP based storage rings and atom beam 
splitters. We prove rigorously that the geometric phase in 
the center of the storage rings proposed in Refs. [H, [l^ 
is always zero. In addition, we find that in the storage 
ring of Rcf. , the spatial fluctuation of the geometric 
phase sensitively depends on the position of the trap cen- 
ter on the "resonance toroid." In the proposals of Refs. 



0, HOi HH, [13 , the fluctuation for the geometric phase 
becomes significantly suppressed when the amplitude 
of the rf-field is large. In the proposals of [H], [l^, the 
fluctuations of the geometric phase also is suppressed if 
the angle rj is set to be k2tt. In the beam splitter real- 
ized with the double wcU potential ARFP [3, [2l|, the 
geometric phase is shown to be zero. 

Our work helps to clarify the working principle of trap- 
ping neutral atoms in an ARFP and the validity condi- 
tions for the various approximations involved. We hope 
our results will shine new light on the proposed inertial 
sensing experiments based on trapped atoms in ARFP. 
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